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Abstract: This paper explores whether Eguchi-Kawai reduction for gauge theories with 
adjoint fermions is valid. The Eguchi-Kawai reduction relates gauge theories in different 
+f numbers of dimensions in the large N limit provided that certain conditions are met. In 

Q-i principle, this relation opens up the possibility of learning about the dynamics of 4D gauge 

theories through techniques only available in lower dimensions. Dimensional reduction can 
be understood as a special case of large N equivalence between theories related by an orbifold 
projection. In this work, we focus on the simplest case of dimensional reduction, relating 

^ a 4D gauge theory to a 3D gauge theory via an orbifold projection. A necessary condition 

t" for the large N equivalence between the 4D and 3D theories to hold is that certain discrete 

^2 symmetries in the two theories must not be broken spontaneously. In pure 4D Yang-Mills 

theory, these symmetries break spontaneously as the size of one of the spacetime dimensions 
shrinks. An analysis of the effect of adjoint fermions on the relevant symmetries of the 4D 
theory shows that the fermions help stabilize the symmetries. We consider the same problem 
from the point of view of the lower dimensional 3D theory and find that, surprisingly, adjoint 
fermions are not generally enough to stabilize the necessary symmetries of the 3D theory. 
In fact, a rich phase diagram arises, with a complicated pattern of symmetry breaking. We 
discuss the possible causes and consequences of this finding. 
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1. Introduction 

In the large N limit, gauge theories can have the remarkable property of volume independence. 
Under some circumstances, one or more of the dimensions of the spacetime in which the large 
N gauge theory lives can be shrunk, while a large set of observables remain unchanged. In 
theories where this volume independence works all the way to zero size, large N volume 
independence opens up the prospect of relating four dimensional gauge theories to lower 
dimensional counterparts that can be studied by techniques only available in D < 4, for 
instance light-cone quantization or quantum mechanical variational methods. Also, it may 
be possible to numerically simulate the lower-dimensional theories at a lower computational 
cost. 

The first version of the observation of volume independence in gauge theories is due 
to Eguchi and Kawaifl], who argued that four dimensional lattice Yang-Mills theory in the 
large N limit is equivalent to a matrix model with no spacetime dimensions (a one plaquete 
model). However, there are conditions that need to be satisfied for the volume independence 
to hold. One of them is that no phase transitions occur as the size of the spacetime volume 
changes. Unfortunately, in the case of pure Yang-Mills, a phase transition does occur as one 
of the dimensions is shrunk below a certain critical size[2]. This is the finite temperature 
phase transition leading to deconfinement, where the center symmetry is broken. Thus, the 
reduction of the spacetime to one plaquete fails in this case. Some attempts were made to 
engineer constructions that would avoid this problem: the quenching procedure[2], the twisted 
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model[3] and, more recently, the double-trace deformation [4, 5]. There are indications that 
the quenching[6] and the twisting [7] procedures fail to protect center symmetry near the 
continuum limit. 

Parallel to these developments, a seemingly different kind of equivalence between large 
N gauge theories was discovered. It was found that one can define 'orbifold projections' 
that relate a 'parent' gauge theory to a 'daughter' gauge theory, which is identical to the 
parent except that all degrees of freedom not invariant under some discrete symmetry of the 
parent theory are left out of the daughter theory. In the large N limit, provided that certain 
conditions are met, there is then a class of 'neutral' observables in the parent and daughter 
theories which are the same in both theories. The parent and daughter theories are then 
termed to be 'orbifold equivalent'. 

Such orbifold equivalences were first discovered in string theory [8], but were quickly 
realized to be purely field-theoretical effects [9]. Proofs of orbifold equivalence at the level 
of perturbation theory can be constructed by analyzing Feynman diagrams[9, 10]. Orbifold 
equivalence and related techniques such as orientifold equivalence have mostly been used in 
the literature so far to relate supersymmetric and non-supersymmetric theories [11]. 

It was pointed out in ref. [13] that large N volume independence can be understood as a 
special case of orbifold equivalence. Volume independence can be seen as a volume reducing 
or a volume expanding transformation that leaves certain observables unaffected, and these 
transformations can be viewed as orbifold projections. In viewing volume independence as 
a special case of orbifold equivalence, one is forced to work with the large N gauge theory 
defined on a lattice. 

For dimensional reduction, one starts with a theory which has one spacetime direction dis- 
cretized on a lattice with T points, and then projects by the discrete translation symmetry Zp. 
Only the translation invariant fields survive the projection, so that the daughter theory is ef- 
fectively defined in one dimension less than the parent. In the volume expansion case[12, 13], 
when one performs 'dimensional reconstruction', one starts with an SU(iVT) gauge theory 
containing some unitary, adjoint scalars, and projects it by a discrete Zr symmetry acting on 
the color indices. If the action of Zr on the different fields is chosen properly, the resulting 
daughter theory is a "moose" or "quiver" SU(iV) r gauge theory with bi- fundamental scalar 
fields. As is well known [14], such theories can be interpreted as theories with one extra 
dimension discretized on a lattice with T sites. The role of the gauge links in the extra di- 
mensions is played by the unitary scalars. The large N equivalence between neutral quantities 
in the parent and daughter theories then amounts to an equivalence between theories "living" 
in different dimensions. 

A non-perturbative proof of large N orbifold equivalence determines the necessary and 
sufficient conditions for its validity [15]. These conditions must be satisfied in both the parent 
and daughter theories, and are i) unbroken center symmetry (confinement) and ii) unbroken 
Zr symmetry, the symmetry used in the orbifold projection. The violation of the first condi- 
tion is what invalidates the original Eguchi-Kawai construction. Eguchi-Kawai reduction is 
valid as long as center symmetry would be unbroken, which is the case as long as the large N 
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theory lives in a volume larger than a certain critical size of order . Numerical results 

support this statement [16]. The proposed modifications of the Eguchi-Kawai construction 
such as the quenched and twisted models, as well as the double-trace deformed models, are 
attempts at preserving the symmetries necessary for large N volume independence (or orbifold 
equivalence) for volumes smaller than this critical size. 

In this paper, we will focus on the simplest case of dimensional reduction, and relate a 
4D gauge theory to a 3D gauge theory by an orbifold projection. We will consider the effect 
of adjoint fermion matter fields on the realization of the discrete symmetries necessary for 
orbifold equivalence to hold so that large N dimensional reduction can work. It has been sug- 
gested that adjoint fermions can prevent center symmetry breaking [13] when a compactified 
spatial direction gets small in a 4D theory 1 . This suggestion was motivated by examining 
the behavior of the effective potential of a traced Wilson loop wrapping a compact spatial 
direction in the AD theory, which is a gauge-invariant order parameter for center symmetry 
breaking. In the limit where the circumference of the compact spatial direction is small com- 
pared to hQQ D , the effective potential can be calculated perturbatively. It was found that 
at one-loop order, the gauge boson contribution to the effective potential favors center sym- 
metry breaking, while the contribution of each flavor of adjoint Majorana fermions is equal 
and opposite to that of the gluons, and favors the preservation of center symmetry. The 
observation of this behavior in a 4D theory motivated the suggestion that 4D YM theories 
with adjoint fermions in the large N limit have volume independence to zero size, so that one 
can construct dimensionally reduced theories, for instance 3D ones, that are equivalent in the 
large N limit to 4D theories due to orbifold equivalence [13]. 

To be sure that large N volume independence really holds to zero size in YM theories 
with adjoint fermions, however, it is necessary to check whether the necessary discrete sym- 
metries are realized appropriately in the dimensionally reduced theory as well as in the 4D 
theory. First, both the 3D and 4D theories should be in a confining phase, which is expected 
to be the case provided that both theories live in large enough spacetime volumes. (Of course, 
as ref. [13] showed, the 4D theory stays in a confining phase even in small spatial volumes 
if adjoint fermions are present.) Next, the discrete Zp translation symmetry along the com- 
pact direction of the 4D theory must not be broken spontaneously. Fortunately, translation 
symmetry is not expected to break spontaneously in gauge theories like the ones that we are 
considering. 2 Finally, it is crucial that the Zp symmetry of the 3D theory must also not be 
broken spontaneously for large N dimensional reduction through orbifold equivalence to work. 
It is this last question that we will focus on in this paper. 

In section 2 we engineer a three dimensional theory that, upon a volume-expanding 

1 Adjoint fermions only help protect spatial volume independence, since they must have periodic boundary 
conditions to protect center symmetry. When fermions have antiperiodic boundary conditions, they do not 
help to protect center symmetry [13]. 

2 At non-zero chemical potential for baryon number in a large N theory, translation symmetry might break 
spontaneously due to the formation of some kind of nuclear matter (for instance something like a Skyrme 
crystal of baryons) at a certain critical value of the chemical potential. 
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orbifold projection, leads to a latticized version of a four dimensional gauge theory with 
adjoint fermions. In section 3 we compute the relevant effective potential and discuss the 
spontaneous breaking of the Zr symmetry used in the orbifold procedure. The minimization 
of this potential is discussed in section 4, and we find that, for most of the parameter space, 
Zr is broken and the large N equivalence fails. We then discuss the origin and consequences 
of this very surprising result. 

2. Dimensional reconstruction and fermions 

In this section, we construct a 3D theory which will be related to a 4D YM theory with adjoint 
fermions, and discuss its relevant discrete symmetries. We then apply a volume-expanding 
orbifold projection to the 3-D theory, and show that the projection produces a 4D theory with 
adjoint fermions. 

2.1 D = 3 theory 

We will start with our 3D theory, which will be the parent from the point of view of the 
orbifold projection. This theory is engineered in such a way as to generate, after a volume- 
expanding orbifold projection, a daughter theory which is 4D, SU(iV) gauge theory coupled 
to one flavor of adjoint Dirac fermions. The AD theory will live on M 3 xS 1 , and the circle will 
be discretized on a lattice with V sites. The fermions will have periodic boundary conditions 
on the circle, so that we are considering a spatial compactification in the 4D theory. 

Since the 4D theory will have fermions and will be defined on a lattice, a naive dis- 
cretization of the fermions will lead to fermion doubling in the continuum limit. To prevent 
doubling, we will engineer the 3D theory to give rise to a 4D theory with Wilson fermions, as 
it turns out to be easiest to work out the orbifold prescription in this case. As is well-known, 
the Wilson term breaks chiral symmetry in the 4D theory, and induces an additive renormal- 
ization of the fermion mass in the 4D theory. This means that we also have to introduce a 
bare quark mass term, so that the 4D theory can have light fermions in the continuum limit 
if the bare quark mass is tuned appropriately against the Wilson term coefficient. 

The 3D theory is a SU(AT) gauge theory coupled to one adjoint unitary scalar and a 
pair of adjoint fermions. We work in Minkowski space with the mostly minus metric. The 
action of the 3D parent theory is 



where a = 0, 1, 2, p a are the 3D Dirac matrices, the D a are covariant derivatives in the adjoint 
representation, x 1S a doublet of two-component 3D adjoint Dirac fermions, r 1 and r 3 are 
the Pauli matrices o"i,o"3 acting on the x flavor space, and (j) £ SU(NT) is a unitary scalar 
field. The constants g, a, T and r will turn out to be, respectively, the AD dimensionless 




(2.1) 
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coupling constant, lattice spacing, the number of lattice points in the fourth direction, and 
the coefficient of the Wilson term. To make the behavior of these fields once we move to 
the 4D theory more transparent, we have written the 3D action above in terms of fields with 
mass dimensions normalized according to the usual 4D conventions. 

The kinetic term of the scalar <j) will generate the gauge-kinetic terms in the X3 direction 
after the orbifold projection. Note that somewhat surprisingly, the would-be Wilson term 
(the term proportional to r/a in the action) has the same mass dimension as the would-be 
X3 kinetic term for the fermions (the term proportional to i/a) in the 3D action Eq. (2.1). It 
will turn out that upon orbifold projection, these two terms will behave differently, and the 
Wilson term will become a dimension 5 operator in the 4D theory, as it must. 

We use the D = 2 + 1 Dirac matrices defined as 



ia 2 , 



(2.2) 



where the g % are the Pauli matrices. 

The action in Eq. (2.1) has a discrete Zp symmetry which acts on the fields as 



A — ► jArf , 



(2.3) 



where uj = e r and 



1 



N 



Ujl 



N 



N 



V 



w r_1 ljv / 



(2.4) 



where ljv is the N x N identity matrix. This Zr symmetry will be the symmetry used to 
define the orbifold projection. 



2.2 Orbifold projection to D = 4 

The volume-expanding orbifold projection amounts to dropping all degrees of freedom in 
Eq. (2.1) not invariant under Eq. (2.3). The only surviving components of the gauge and 
fermion fields are in the N x N diagonal blocks, while the surviving degrees of freedom of the 
scalar field are in the N x N one-off-diagonal blocks. 
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X2 
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A, 



X 



X3 



(2.5) 



V 



A r / 



V 



Xr / 



/ 



The action of the daughter theory is, up to a factor, the orbifold-projected action of the 
parent (3D) theory S3: 



Ss[A,fax] 



1 r_1 

S D [A 1} ■ ■ ■ A r , fa • • • r , Xl • • • Xr ] = - £ Sp[ 7 n A 7 nt , y^y* ^y^t] 



n=0 



(2.6) 

When valid, the orbifold equivalence states that, to leading order in 1 /N, the correlators of 
neutral operators in the parent (that is, operators invariant under the the Zr transformation 
in Eq. (2.3)) agree with the correlators in the daughter theory. An example of a neutral 
operator in the 3D theory is tr^> r , which is a Wilson loop wrapping the compact S 1 direction 
from the point of view of the 4D theory. 

In the case of the action Sp in Eq. (2.1) we have the daughter theory action 



r r 

S D = d 3 x aJ^Tr __Lf^F^ + ^\D a fa\ 2 + X i(ip a D a + mr 1 )^ (2.7) 
i=i L 9 a 9 

~^Xi4>\ - <t>lxi+i)T"\xi(t>i + 4>iXi+i) + 7fa(xi4>t - ^Ixi+iWix^i - faxi+i) 



where 

D a fa = d a <j>i + i{A l a fa - faA^ 1 ), 

D aX i = d a Xi + i[Al,Xi]. (2.8) 

The action in Eq. (2.7) is the action for a 4D gauge theory discretized in the compact x 3 
direction, with the unitary scalar fa as the i-th link variable in the x 3 direction. In fact, 
writing fa = e MA3 W we see that the kinetic term for is a discretized version of F a 3: 

D a fa ~ d a fa + iA a fa - ifaA 1 ^ 1 

~ m fd Q A 3 - AL+1 ~ AL + i(A> 3 - A'A^ 1 )) ~ iaF a3 , (2.9) 

Similarly 

r 

Y,(Xi+ir 3 4xifa - XiT 3 faXi+i4) - 2ia X iT 3 D 3X . (2.10) 
i=i 



- 6 - 



and 

r 

^(Xi+iTVPx^ + XiT^Xi+lti ~ 1XiT X Xi) ^ -a 2 X T\Dz) 2 X . (2.11) 

i=l 

We now define the four-component AD Dirac spinor ip by 

i- = ( , (2.12) 

and the AD Dirac matrices in the chiral basis in terms of the 3D Dirac matrices as 

with = (1, (?) and cr^ = (1, — a). Essentially, the 3-D Dirac spinors become AD Weyl spinors 
in the orbifold construction. 

With the identifications above, we see that in the small a limit the action in Eq. (2.7) 
becomes 



S d = / d x tr 



._L F 2 + ^(yfDp + m)i/> - ^D^Dsi/j + 0{a 2 ) 



(2.14) 



As promised, the orbifold projection takes a theory defined in 31? and generates a AD 
theory, with the Wilson term becoming an irrelevant dimension 5 operator, as it must. The 
AD theory has gauge group SU(iV), in contrast to the SU(iVT) gauge group of the 3D parent 
theory. The spatial extent of the reconstructed x 3 direction equals L = Ta and thus can be 
adjusted by tuning the parameters of the 3D theory. 

It is possible to adapt the construction above to accommodate AD theories with Majorana 
fermions, such as N = 1 Super Yang-Mills theory, which has one flavor of massless Majorana 
adjoint fermions 3 . To do this we can impose the condition cr 2 (x 1 )* = °" 3 X 2 i n Eq. (2.12), so 
that the 4D Dirac spinor written in terms of 3D spinors has the form 

*- («£.)■ (2 ' 15) 

and becomes a AD Majorana spinor. Setting c 2 (x 1 )* = °" 3 X 2 ; the action S3 becomes 



^Majorana = ^ j ^ tr 



— F 2 + -^\D a cP\ 2 + 2xip a D aX + m(xW + X T <r 2 x)(2.16) 
2g z 2g z a z 



-{bcA ] ]{x\4>} - ix*,^]{x*,<P}) - ^{[x*A ] ][x\4>] + [x,Hx,<t>]) 

a 2a 



where now x does not carry any flavor indices. 

J Supersymmetry would be broken by the lattice discretization. 
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3. Large N equivalence and the effective potential 



For the 3D theory to be orbifold equivalent to the AD theory in the large N limit, two basic 
conditions must be met. First, both theories must be in their confined phases, and second, 
the symmetry that defines the orbifold projection, in our case the Zp symmetry in Eq. (2.3), 
must remain unbroken in the parent theory. We expect the 3D theory to be in a confining 
phase as long as it is defined in a large enough 3D volume. The AD theory is expected to be 
in a confined phase when L = Ta is large enough, and also when L < Aqq D because of the 
presence of the adjoint fermions[13]. The crucial question that will determine whether large 
N orbifold equivalence between the two theories is valid is whether or not the Zr symmetry 
of the 3D is spontaneously broken. 

A set of gauge invariant order parameters for the Zp symmetry is comprised of the expec- 
tation values of the traces of powers of (f>, (tic<f> k ), for k = 1, • • • , T. The information encoded in 
the r quantities (cj) k ) is also encoded in the eigenvalues of (p (but not in their ordering, which 
is a gauge-dependent quantity). We will investigate whether Zp is spontaneously broken by 
calculating the effective potential for the eigenvalues of 0. Actually, we will do a more general 
calculation. We can diagonalize <f> through a gauge transformation Q as (j) = QcftoQ,^ , where 
00 is diagonal: 



/ e Hpi 



\ 



y e ltpNr ) 

In order to test for Zr breaking, it would be enough to consider a 4>q with only T distinct 
eigenvalues, so that <f>o has a 'block diagonal' form if they are properly ordered. We will, 
however, consider the more general case where all eigenvalues can be different, and test for 
the breaking of the full Z^vr symmetry defined by Eq. (2.3) with uj = e 2m / Nr anc [ 



(3.1) 



/l 



\ 



NT-1 



(3.2) 



V w / 

The effective potential for the eigenvalues of 4> will then determine how the Z^vr symmetry 
and its Zr subgroup are realized. Of course, the computation of the effective potential for 
{trcj) k ), or equivalently for the distribution of the eigenvalues of (f>, is in general a difficult 
non-perturbative problem. However, we expect the relevant scale for the renormalization of g 
to be the size L of the "extra" dimension 4 , as this is the case for gauge theories in M 3 x S . For 
L = aT <C Aq^ d , where Aqcd is the strong scale of the AD theory, we expect a semiclassical 



The only other scale available in the 3D theory is a < L. 
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approximation to be valid, and we will compute the one-loop effective potential for the "Lnv 
symmetry. 

The one-loop effective potential will be given by the small (quadratic) fluctuations around 
4> = 4>q. A perturbative calculation like this requires gauge fixing, and we will use gauge. 
We parametrize the (j) field by <fi = e tG 4>oe~ tG , where the G's are the would-be Goldstone 
bosons. 

Gauge-fixing amounts to the addition to the lagrangian of 

C^C-\FF^ + c( 5 -^)c, (3.3) 
where c, c are the Faddeev-Popov ghosts. We choose the gauge fixing function F to be 

F=^ 7 =[dA--^GM (3.4) 
After gauge fixing, the part of the lagrangian quadratic in A, dcpo, G and bilinears in c, c, 

is 

^Cd™ t!C = "^2 tr ^2d a Apd a A^ - 2(d.A) 2 + 1[A, to][A,^] - ±[dA,<M[dA,</>l]^ 

+trc(-a 2 - 4)[c, M + ^-^tria^ol 2 

+tix(ip a D a + mr l ) X - ^[ X , 4}r 3 { X , M - ^ [ X , ^[x, <h] • (3-5) 

The would-be Goldstone bosons in G pick up a gauge dependent mass tvlq ~ £/a 2 , the ghosts 
pick up a mass m 2 ~ £/a 2 and the longitudinal components of the vl that gets higgsed also 
pick up a mass m\ ~ £/a 2 . The transverse components of the higgsed A's (we shall call them 
W's) pick up a real mass mw ~ 1/a and the fermions x a ^ so get a mass m x ~ m + 1/a. The 
precise values of these masses depend on the value of 4>q- To find them, let us parametrize 
the adjoint Hermitian fields as A = A a X a , with trA a A b = ^S ab , and similarly for G, x, c and c. 
We need to evaluate terms such as tr(A a 0QA fe (/>o)> which is easiest to do by considering the 
decomposition of SU(NT) into its SU{2) subgroups. It is not hard to show that if A a ,At 
belong to different SU{2) subgroups, or to the Cartan subalgebra, they do not make a 4>o 
dependent contribution to the effective potential. This means that we only need to consider 
a sum over SU{2) subgroups with A a , Aj, belonging to the same subgroup. 

To show how this works, we will temporarily work with NT = 2; the generalization to a 
general SU(NT) gauge group is straightforward and will be done at the end of the calculation. 
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We will need the following traces: 
tr(A°4A%) = 

tr([A a ,4]{A 6 ,<M) = tr(A a 4AVo - A> A 




sin(^i - <fo) (T 
sin(0i - <j> 2 ) 

0, 



tr([A a ,4][A b ,0 o ]) = tr(A>T A ^ + A> A b ^ - 2A a A b ) 




(3.6) 

1 o o\ 

cos(0i — (j) 2 ) — 1 . 

0y 



Using Eq. (3.6) we see that only the 1 and 2 components of the fields are affected by 4>q 
and pick up a 4>o dependent mass, so we will drop the 3 rd component from now on. We arrive 
at 



p fixed 

j-i ■'-'quadratic 



- A (daA^A^ - (d.A 1 ) 2 + i(cos(^! - ip 2 ) - I) A 1 A' 

9 i=l,2 ^ a 



(cos(<^i - v? 2 ) - i)9.^a.^ 



■^(cos(^ - - 1) £ (sew - lcf*c*) 



+(cos(^ 1 - ¥ > 2 )-l) ]T c l (-<9 2 - 4^ 



i=l,2 



i=l,2 



J] x\ip a D a + mr 1 )^ - ^ sinfai - ^(xVV ~ X 2 tV 



-1^3, 2 -2J li 



+ — (cos(^i - y> 2 ) - l)xVY- 
2a 



(3.7) 



where the superscripts in x 1,2 denote adjoint color, not flavor. 

The propagator for the color components 1 and 2 of the gauge fields can now be shown 
to be 



1 



, 2 l-cos(y 1 -y 2 ) 

ru — 



+ 



i 



k 2 (i-cos^x-^))^ 2 -^)' 



(3.8) 



The two physical, transverse degrees of freedom have a ^-independent mass. The unphysical, 
longitudinal component has a ^-dependent mass and wave function normalization. This mass 
(and wave function normalization) is identical to that of the would-be Goldstones and of 
the ghosts. Counting the unphysical degrees of freedom (d.o.f.) we have 2 gauge d.o.f., 2 
Goldstone d.o.f., and 2 ghost d.o.f. As the ghosts contribute to the effective potential with an 
additional factor of —2, the contribution of the unphysical degrees of freedom to the effective 
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potential cancels out. We are left with the physical degrees of freedom: 4 gauge field d.o.f. 
(from two adjoint colors and two polarizations) and 4 fermion d.o.f. (from two adjoint colors 
and two "flavors"). The Higgs field contribution does not depend on cj>Q, so we can ignore it. 
The one-loop effective potential will depend on cj>o only through the value of these physical 
masses and will, consequently, be £— independent. 

The bosonic contribution to the effective potential can now be shown to be 

3/2 



' lH ' log ( lr 2 , sin 2 ' ~ - 2 



■ 2 ( <P1 ~ ¥2 

2 sm 1 



(2tt) 3 to V « 2 V 2 J J 6vra 3 
where we have used the dimensionally regulated and minimally subtracted integral 

± I lo g (k 2 - M 2 ) = T ±- (M 2 ) 3 / 2 . (3.9) 

At this point, one might worry that the integral above might have an infrared divergence for 
some distributions of eigenvalues, since the M 2 of the bosons vanishes for (p% = ip2- This 
is a serious issue in a full Eguchi-Kawai reduction to a D = one-plaquette model, where 
the contribution of the gauge bosons to the one-loop effective potential becomes singular, 
and one must examine higher order corrections to the effective potential to make reliable 
conclusions [17] 5 . However, the integral above for D = 1,2,3 is IR- finite and well-behaved as 
M — > 0, so that this issue does not arise for D > 0. 

The dispersion relation for the fermions can be found in a similar way, where we continue 
to work with the iVT = 2 example for clarity. We find 

1 x f d*k 1 /# + ( m + £(1 -008))^ -^sinr 3 \ |V\ 

Ta quadratic J ( 27r )3 2 U X '\ L s i n r 3 # + (m + ^(1 - COs))^ M X 2 



M 

(3.10) 

where we define the shorthand notation cos = cos(<^i — (^2) and sin = sin(< / c 1 — ^2) and the 
indices on x 1 ' 2 refer to color space. The fermion matrix can be written as 

M= (ft+ (m+ -(1- cos)r 1 )I + -sinr 3 S 2 (3.11) 
a a 

with I and S 2 = a 2 acting on color space. The fermion contribution to the effective potential 
is given by the Tr Log(M) and can be calculated as 

d 3 k 



/d Ik 
j^-^tT co i or tiyi avor tr spin log M 

= 4 |^03 log ^ 2 -((m+^(l-cos)) 2 + (^sin) 2 )) 

-1 3/2 



37ra 3 



[ma + 2rsin 2 (^— -^)) z + sm 2 (ipi - (p 2 ) 



(3.12) 



3 We thank Larry Yaffe for bringing this issue to our attention. 
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Above we used relations like 



tr spin log(^ + A) = \og{k 2 - A 2 ). (3.13) 
Adding up the gauge boson and fermion contributions to the effective potential we find 



V. 



SU(2). 



1 



eff W) 



67ra 3 



4V8| sin( 



3/2 

11 J (ma + 2r sin ( )) + sm (^1—992) 



It is straightforward to generalize the calculation outlined above to the SU(iVT) case, 
with the result that 



i<j=l,...,NV L 



V2\ sin( 



2 

This is the main result of this section. 



_ 3 / 2 " 

(ma + 2rsin 2 ( ^ - — )) 2 + sin 2 (c/?j - ipj) 

(3.14) 



4. The phase diagram for Z Nr symmetry 

A few observations about the effective potential in Eq. (3.14) can help us gain insight into 
the result. First, let us consider the purely bosonic case. In this case the effective potential 
reduces to 

-Is E i-(^)i 3 . ("J 

i<j=l,...,NT 

It is clear that the effective potential is minimized when all of the eigenvalues are degenerate 
and (tr <^>q) 7^ 0. This ground state breaks Zjvr maximally, down to nothing. In particular, 
the Zr symmetry used to define the orbifold projection is spontaneously broken, which implies 
that the equivalence between the large N dynamics of the 3D and AD theories is lost. This 
conclusion agrees with the standard belief that the Eguchi-Kawai reduction for pure YM 
breaks down if one (or more) dimensions are shrunk down to one plaquete. 

Let us now look at the fermion contribution to the effective potential. First, note the 
surprising fact that the Wilson term coefficient shows up in the effective potential, and does 
not disappear in the small a limit. This can be traced back to the fact that the 3D progenitor 
of the 4D Wilson term shows up at the same order as everything else in the 3D action, even 
though the Wilson term is an irrelevant operator in the 4D theory. The effective potential of 
the 3D theory thus depends strongly on the value of the Wilson term coefficient. 

Next, note that the fermion mass will be additively renormalized from its bare value by 
loop contributions due to the couplings of the fermions to the gauge field. Those effects are 
higher order in g 2 and are not reflected in Eq. (3.14). The renormalization of the mass would 
show up at two-loop order in the calculation of the effective potential. Thus, at the order to 
which we are working, the parameter m in Eq. (3.14) can be viewed as the "fermion mass" 
(with the obvious issues relating to the proper definition of fermion mass in a confining theory 
applying). 
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Now we must ask what region in the parameter space of the 3-D theory can correspond 
to a continuum AD theory, assuming that the orbifold equivalence is valid. The continuum 
limit in the AD theory corresponds to sending a — > 0, V — > oo in such a way that Ta = L 
is fixed. This means that the lattice spacing a must be much smaller than any other length 
scale in the problem, and in particular: 



This means that we should only expect to have our 3D theory to be equivalent to a continuum 
AD theory with fermions in a region in parameter space where ma <C 1. 

The general effect of the fermion contribution to V e fi naively seems to counteract the 
symmetry breaking effect of the gauge fields, as the bosonic and fermionic contributions have 
opposite signs. For instance, for the special case ma = and r = 1, where we might hope that 
the equivalence should work, the potential between any two eigenvalues is minimized when 
they are maximally separated {<pi — <pj = ±7r). That does not mean, however, that the Zjvr 
symmetric distribution of eigenvalues (where e l ^ n = e t2nm / Nr ( n = 1, • • • ,NT) is favored, as 
the symmetric configuration does not minimize V e s for each pair of eigenvalues. V e g can, 
instead, be minimized by having half of the eigenvalues equal 1 and the other half equal to 
— 1. The energy lost by having pairs of identical eigenvalues is more than compensated for 
by having even more pairs with the maximal separation. As a result, the Zp symmetry used 
in the orbifold projection breaks spontaneously, and we see that surprisingly large N orbifold 
equivalence is not valid for ma = 0, r = 1, and the 3-D theory is not equivalent to a continuum 
AD theory. 

For generic values of ma and r, the minimization of V e g is a complex problem and better 
dealt with numerically. The result of such a minimization is summarized in Fig. 1. There we 
can see a rich phase diagram as a function of ma and r, particularly in the region away from 
1 3> ma, where we are probing lattice-scale physics in the AD theory. The purple region is the 
Zjvr symmetric phase; other colors denote less symmetric phases. Note that at small values 
of ma, where we might hope the orbifold equivalence might hold and our 3D theory would 
correspond to a continuum 4D theory with fermions, the fermion contribution is unable to 
prevent the breaking of Zr and the large N equivalence between the 3-D and the AD theory is 
lost. Even at ma = and r = 0, where the AD theory would have Dirac fermions with lattice 
doublers, which should help stabilize the center, the eigenvalues are in the totally collapsed 
phase, and Zr is broken. A similar behavior was observed in one-loop calculations of effective 
potentials in matrix models of M-theory [18]; the breaking of the symmetry there, like here, 
was interpreted as a failure of the model to "generate" spacetime dimensions. 

5. Discussion 

The breaking of Zr may seem at odds with similar calculations done in the past [19], for 
instance the calculations suggesting that adjoint fermions protect orbifold equivalence and 
allow large N dimensional reduction. However, there is no contradiction: we would only 



ma <C 1 . 




- 13 - 




-l.S -I.U -0.5 0X1 0.5 1.0 1.5 



Figure 1: Phase plot of one-loop effective potential with axes r vs. am in the large AT limit. 
Configurations of the AT eigenvalues distributed at 1 , . . . , i . . . , AT points evenly spaced around the 
circle were tested against each other, and the configuration with the lowest energy was determined 
at each point in the plot above. We refer to the configuration with the eigenvalues distributed on i 
points on the circle as the Zi configuration in the plot legend above. The purple region is the region 
where the full Z^v symmetry is unbroken. 

expect the 3D effective potential to match the 4D effective potential for the trace of a Polyakov 
loop wrapping the compact direction if the two theories were large N equivalent. If the two 
theories are not equivalent, the 3D and 4D effective potentials do not have to agree even 
if they are effective potentials for the same order parameter. As it happens, however, the 
effective potential that we have calculated in the 3D is for different order parameters than in 
the 4D calculations. 

In the 4D calculations in the literature, SU(A r ) Yang-Mills theory [13, 19, 20], possibly 
with some matter content, was taken to be defined on IR 3 x S 1 , with the size of the compactified 
dimension equal to L, and the calculations were done in the continuum limit. To determine 
whether large A?" volume-independence holds down to volumes small compared to the strong 
scale, the fate of the Z^r center symmetry was determined by looking at the effective potential 
for the Polyakov loop wrapping the compact spatial direction. These calculations have shown 
that adjoint fermions protect the center symmetry at small volumes. 

Such 4D continuum calculations do not however answer the question of whether in a large 
N lattice theory with adjoint fermions (discretized along at least the compact direction), one 
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can discard all but one lattice site along one direction and obtain a large N equivalent theory, 
which is the relevant question for determining whether orbifold equivalence holds and large 
N dimensional reduction works. For that, one must examine the realization of the relevant 
discrete symmetries in the higher and lower dimensional theories, and it is precisely the latter 
question that our calculation of the effective potential in the 3D theory addresses. Our results 
suggest that the relevant symmetries in the 3D theory break spontaneously, so that one does 
not obtain a large N equivalent theory in a full dimensional reduction even in the presence of 
light adjoint fermions. 

However, there remains a major puzzle. There is a region in parameter space, which is 
colored purple in the phase diagram in Fig. (1), where the full Z^r symmetry is unbroken, 
and orbifold equivalence appears to hold, at least at first glance: all of the conditions for 
large N orbifold equivalence appear to be met in this region of parameter space. However, in 
this region the value of the fermion bare mass (equal at this order to the renormalized mass) 
is of the order of the cutoff ma ~ 1 and the fermions should effectively decouple from the 
low energy dynamics of the 4D theory. The 4D theory should then be pure Yang-Mills in the 
continuum limit. That is, again, contrary to any expectation. 

To get further insight into our results, we can compare a quantity that can be calculated 
in both the 3D and 4D theories, which is the effective potential for a Wilson loop wrapping 
the compact direction in the 4D theory. For small enough L <C 1/A, when it is legitimate to 
compute the effective potential in the semiclassical expansion, the effective potential can be 
calculated in both the 3D and 4D theories. 

In the pure Yang-Mills 4D theory, the effective potential for the Polyakov loop can be 
calculated by summing over the small fluctuations of the gauge field around a background 
gauge field with a non-trivial Polyakov loop: 



A, t = A 3 6l+A„ A 3 = diag(0i/L, • • • ,0 N /L), fi 3 = P e< J" <fa8 * = 



e i8 2 



y e i0N J 

(5.1) 



The result, after gauge-fixing similar to Eq. (3.14) in 3D, is [13, 19, 20] 



lAA /■,„,, /27m + ^-0A 2 \ 1 N °° ' 



j 



«»~z££y PT^J -fi££ 

i,j=l n=l J \ V 7 / i,j=l 71=1 V 

1 o \ e i-°j\\ 1 cos(n(0j - 6j)) 

z,j=l i,j=l 71=1 

oo 

■^E^l^sf, (5-2) 

71=1 

where we used the ("-function regularization to define the sum over the Kaluza-Klein modes, 
and used Hurwitz's formula in the next to last equality. The result can be written in two 
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useful ways [13] , one making the sum over the winding numbers of the Polyakov loop explicit 
and the other making the sum over the eigenvalues explicit: 

n=l 

v*d = -— ^ (^f - x> - o 3 m - o 3 ] - 2,r j . (5.3) 

where [a] = a mod 2ir. It is clear from these expression that the potential is minimized when 
all eigenvalues collapse to one common value. Each flavor of adjoint Majorana fermions makes 
an equal and opposite contribution to that of the gauge bosons to the effective potential, so 
that adding one flavor of Dirac fermions, or equivalently two flavors of Majorana fermions, 
will simply have the effect of flipping the sign of the effective potential above, making it 
repulsive. The resulting potential favors Zjy symmetric distributions of eigenvalues of the 
Polyakov loop. 

We can also calculate the effective potential for the Polyakov loop from the 3D theory. 
In fact, the calculation of the effective potential we have done above is general enough that 
it already gives the answer in the appropriate sector. If we consider a <f>o that takes the 
block-diagonal form 

/e^l r \ 
00= ••. , (5.4) 

\ e iVN lv) 

where lr is the r x V identity matrix, Tifi = ft, then r = Q3 and, in the region of parameter 
space where the large N equivalence holds, the effective potentials for the eigenvalues of 12 3 
calculated in the 4D and 3D theories should be the same. Let us look at the effective potential 
for tr as calculated from the 3D theory by making making the appropriate substitutions 
in Eq. (3.14). In the continuum limit, where T — > 00, a — > 0, with Ta = L fixed, Eq. (3.14) 
reduces to 

1 N 

V eS (fl 3 ) = £ (V2|ft - 9j\ 3 - (8(mL) 3 + 12mL(amr - l)|ft - 6 3 \ 2 )) . (5.5) 

i<j 

Note that the term with the Wilson term coefficient disappears in the continuum limit of 
a — ► 0, as it must. However, this effective potential does not match Eq. (5.3) for any value of 
m. In particular, the two results clearly do not match for m ~ 0.7 /a, r ~ 1, which is in the 
region in parameter space where Zr symmetry is preserved, and large N orbifold equivalence 
would naively be expected to hold. Probably the most important difference between the two 
effective potentials is the lack of anything like a sum over Kaluza-Klein modes in the 3D 
theory, while such a sum appears naturally in the AD calculation of the effective potential. 

In our calculation of the effective potential of the 3D theory, we took the scalar field (j) 
to be close to a block diagonal form. But it was the one-off-diagonal terms in Eq. (2.5) that 
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generated the derivative terms in the "extra" dimension. This means that we had essentially 
no chance to see derivatives of the KK modes arise in our computation, and thus each of the 
terms in the sum over the KK modes is independent of the KK mode number n. The sum 
over n reduces in our case to a factor of T, with the \6i — #j| 3 dependence agreeing between 
the 4D and the 3D calculation. Of course, if large N orbifold equivalence were to hold, the 
behavior of gauge- invariant observables in the common sector of the two theories should agree 
regardless. The fact that the effective potential for center symmetry breaking calculated in 
the 4D theory and in the 3D theory do not agree even in the region of parameter space where 
Zr symmetry is unbroken suggests that the equivalence is breaking down for some other 
reason that is not currently understood. 

One way to reconcile the continuum 4D calculations with our construction is the follow- 
ing 6 . The 4D calculations suggest that L can be reduced arbitrarily as long as L > o, that 
is, within the region of validity of a continuum approximation. A theory like the 3D theory 
described here can be seen as the result of reducing L all the way to L = a (even though what 
we consider in this paper is dimensional expansion from 3 to 4 dimensions and L = Ta > a 
at all times). From that point of view, the results stemming from our 3D theory are lattice 
artifacts. That would suggest that a YM theory with adjoint fermions may have large N 
volume independence which works for arbitrarily small L in physical units, but the resulting 
volume-reduced theory must still have a large number of lattice points in the compactified 
direction, which look like a large number of independent fields form the 3D point of view. 
Viewing such a "reduced" theory as being 3D does not seem to provide any advantage over 
just working with the original 4D theory. 

The challenge of unraveling this situation is left for future work. A working orbifold equiv- 
alence allowing large N equivalence between interesting theories such as Yang-Mills theory 
(with matter or without) in different numbers of spacetime dimensions would be extremely 
useful. Unfortunately, as we have seen in this paper, it is not yet clear how to make such an 
equivalence work. 
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6. Postscript 



After the completion of the present work and partially triggered by it, Bringoltz[21] analyzed 
similar questions from a different perspective. His analysis sheds some light on the issues we 
discuss above, and we would like to comment on it here. 

Ref. [21] pointed out that the 3D reduced model in Eq. (2.1) is non-renormalizable in 
perturbation theory, and that the linear divergences that arise in calculating the one-loop 
effective potential for Zjvr breaking 7 must be absorbed by counterterms coming from certain 
relevant operators that should be added to the action in Eq. (2.1) and not, as we erroneously 
assumed before, by terms already contained in the action. For instance, in the bosonic sector 
the effective potential given in Eq. (3.9) evaluated with a ultraviolet cutoff A is proportional 
to 

^3 j\kk*Io g (k>-M>) = -^- i ^ + 0(A-i), (6.1) 
where M 2 = \ sin 2 ( y ' 2 <fj J , and we have discarded a field-independent divergent constant. 



The linear divergence is absorbed by 

di / d 6 x\tv^ (6.2) 



and not, as previously assumed, by a wave function renormalization. When fermions are 
included two counterterms are required at one-loop 



di J d 6 x\tr(j)\ 2 + d 2 J drx\tc<f>T ■ (6-3) 

It seems then that an interesting question is whether the theory in Eq. (2.1) supplemented 
by these two terms can be equivalent to a four-dimensional gauge theory. For this to be true 
two conditions need to be met. First, Zjvr must not break. Secondly, the theory should be 
orbifold projected into a theory equivalent to the four-dimensional gauge theory of interest. 

The coefficients d\, d 2 contain regulator-dependent divergent pieces, which serve to cancel 
the one-loop UV divergences, and also contain arbitrary finite pieces that contribute to the 
effective potential at tree level. The values of the finite pieces can affect the realization of Zjvr 
symmetry. Thus, to check the first condition above, we added the two double trace terms in 
Eq. (6.3) and varied their coefficients while checking the center symmetry realization using 
the methods used above in Sec. 4. We were unable^ to find double-trace term coefficients 
that would protect center symmetry in the theory with Nf ^ 0. As is illustrated in Fig. 2 
for Nf = 1, the center symmetry breaks to Z3, even when the UV-finite parts of the double 
trace coefficients are large 9 . 

This result is expected by general considerations. The stabilization of center symmetry 
using double trace terms was discussed extensively in Ref. [4] . As argued there, it is expected 

7 Refs. [21, 22] note that if the "extra" dimension is shrunk to two or more points (instead of one), the 
linear divergence we discussed here disappears. This indicates that a continuum 3D model with two unitary 
adjoint scalars should not have the linear divergence we are discussing. This is essentially an example of the 
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Figure 2: The figure on the left shows the Z^r phase diagram following from numerical minimization 
of the one-loop effective potential with only the first two double-trace terms included with coefficients 
e?i, g?2- The figure on the right shows the phase diagram with the full tower of double trace terms, with 
all double trace terms having the same coefficient d, with varying bare fcrmion mass m and Wilson 
term coefficient r = 1. In both figures NT = 25 and Nf = 1. 



that if the coefficients of the terms in Eq. (6.3) are chosen appropriately, one can prevent 
center symmetry breaking with titfi and tr0 2 as order parameters. However, no matter what 
their coefficients are, these two double trace terms would not generally be expected to prevent 
the breaking of center symmetry with higher powers of </> as order parameters, for instance 
breaking with tr</> 3 as an order parameter. (This is breaking pattern referred to as Z3 in 
Sec. 4). As argued by Unsal and Yaffe[4], to fully protect center symmetry using double trace 
deformations, one would need to include a whole tower of double trace terms, with all powers 
of 4> up to AT/2: 

AT/2 

d k I d z x\t^\ 2k . (6.4) 
fe=i J 

This suggests the addition of the full tower of double-trace operators to the action, whether or 
not their presence is required by renormalization at one-loop. Once the full tower of double 

Little Higgs construction [23] popular in electroweak symmetry breaking model building. 

8 The exception is if NT is prime and one allows the double-trace term coefficients to scale with NT, which 
is not what one wants to do for large N equivalence. 

9 The situation for Nf = is somewhat different at one loop, since a direct numericai check shows that 
the two double trace terms in Eq. (6.3) can indeed prevent center symmetry breaking. However, we do not 
expect this to hold beyond the one-loop approximation, and the full tower of double trace terms is likely to be 
necessary to protect center symmetry beyond one loop even for Nf = 0. 



- 19 - 



trace terms in Eq. (6.4) is included, center symmetry can easily be protected by choosing 
appropriate coefficients for the double-trace terms, as is illustrated in Fig. 2 10 

We now turn to addressing the second condition: whether the theory with the added 
double-trace operators orbifolds into a theory equivalent to a four-dimensional gauge theory. 
It can be argued that the theory with the added double-trace terms and dk ~ 1/L 3 = l/(Ta) s 
is equivalent to a theory without them on two accounts in the large volume limit. First, in 
the large L limit the deformation vanishes. Second, as argued Ref. [4], their contribution to 
the dynamics at large volume in the large N limit are suppressed by 1/N 2 . In our case, the 
double-trace terms orbifold into 



where trfi = tr^> r is a traced Wilson loop wrapping the compact direction in the 4D theory. 
Recall that terms like Itr^l 2 ,^ ^ {mT,m S N} do not survive the volume-expanding orb- 
ifold projection. Thus it is clear that the two double trace terms whose coefficients contain 
regulator-dependent UV divergent pieces at one loop do not appear in the orbifolded theory. 
Also, it is clear that at one loop, the coefficients dk with k > 2, some of which survive the 
orbifold, are independent of the regulator of the 3D theory. Since we want dk to make contri- 
butions to the one-loop effective potential of comparable size to the other terms in the action 
(that is, finite in the large T and large N limits), we must generically take dk = A/a 3 , where 
fk are numbers independent of T and N. 

The resulting four- dimensional double-trace terms clearly do not vanish in the large 
L = Ta limit, and the first argument given above for the irrelevance of the double-trace terms 
in the large- volume 4D theory no longer applies. However, the second argument, involving 
the 1/N 2 suppression, will still work provided that the large N limit is taken before the 
continuum limit and before the large volume limit: the limits do not commute. With this 
ordering of limits, once the values of coefficients dk are chosen appropriately, the 3D theory 
will be orbifold-equivalent at large N to 4D YM with adjoint fermions with double trace 
deformations, and at large volume L and large N, this latter theory will be equivalent to the 
desired 4D theory, which is simply 4D YM with adjoint fermions. 

The upshot of the above discussion is that to prevent center symmetry breaking and have 
a 3D theory equivalent to a 4D gauge theory in the large N limit it is necessary to add the 
full tower of double-trace operators to the three-dimensional action and rely on the 1/N 2 
suppression of the double-trace operators in the large N limit. 

Ref. [21] approached the problem not from a continuum theory as we did in this paper 
but worked with a lattice version 11 of Yang-Mills with adjoint Wilson fermions with a variable 

10 Interestingly, one can also preserve center symmetry by putting the full tower of double trace terms except 
for all terms that contain 4> kT , where k are integers. This suggests that one can protect the center symmetry 
with a double trace deformation that does not survive the orbifold projection. The implications of this curious 
result are left to future work. 

n See also Ref. [24] for a related numerical study of a reduction of the 4D theory to one plaquette. 




(6.5) 
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number of sites. When the number of sites in all directions is large, taking the continuum limit 
one obtains continuum 4D Yang-Mills theory with adjoint fermions, plus higher dimensional 
operators. On the other hand, when the number of sites in the three non-compact directions 
is large, but there is only one link in the compact direction, the continuum limit gives the 
theory discussed in the present paper, with the link in the compact direction denoted by </>, 
supplemented by the double-trace operators in Eq. (6.3). In the construction of Ref. [21] 
there are not explicit double-trace operators in the action. However, if one desires to discuss 
the continuum limit of the lattice theory, one can view our continuum theory as the low 
energy effective theory for it. From this perspective, the numerical values of d±,d2 in the 
effective theory are determined from the lattice theory, and are dependent on the particular 
lattice regularization used. Working in perturbation theory it was shown that two double- 
trace operators are indeed generated at one-loop level. As argued before, they are not enough 
to suppress Zjvr breaking. 

At first glance, it might appear that there is some tension between these results and 
that of Ref. [21]. In that work, Z^vr breaking was tested in the lattice-regularized theory 
with adjoint fermions (which assigns some particular values to the double trace coefficients in 
Eq. (6.3)), and it was argued that Z^vr did not break at one loop. This conclusion was based 
on comparing the potential energy of the Zi,Z2, and Zjvr symmetry realizations, and the 
Zat symmetric realization was found to be preferred to the other two tested. If, as suggested 
by our results above, the center symmetry in the lattice-regulated theory breaks with tr</> 3 
with an order parameter, this would not have shown up in the analysis of Ref. [21]. If, on the 
other hand, Z^rr breaking does not occur after all in the lattice-regulated theory, consistency 
with our results above would suggest that the lattice-regulated theory with fermions actually 
generates the whole double-trace tower at one loop in the continuum limit, rather than just 
the first two terms considered in Ref. [21] . 

We now finally return to the question raised in the title of this paper: do fermions save 
large N dimensional reduction? Working with the 3D continuum theory, which is the choice we 
investigated here, the answer is no. The preservation of center symmetry depends entirely on 
the coefficients of the double-trace terms in the action, which generically have UV-divergent 
pieces that depend on the choice of regulator, and finite parts that can be chosen freely by 
hand from the EFT point of view. We see that while it appears that (happily) large N 
dimensional reduction can be saved, fermions do not participate in the rescue: The fermions 
themselves are not enough to protect center symmetry of the 3D theory, and just turning on 
adjoint fermions is not enough to make the 3D and 4D theories orbifold equivalent. 

The answer to the same question is more subtle in the lattice version of the 3D theory 
which was investigated in Ref. [21]. There, it was argued that the lattice regulator of the 
theory with Wilson fermions effectively assigns certain UV-divergent values to the coefficients 
of the first two double-trace terms that turn out to be sufficient to protect center symmetry 
in the 3D lattice theory with one site in the "extra" direction. If this is indeed the case, 
the center symmetry of the 3D lattice theory is saved by UV effects associated with the 
fermion discretization, rather than the infrared physics of the fermions, as might have been 
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expected given the results found in the analogous 4D continuum theory in Ref. [13] with 
adjoint fermions in small volumes. 
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